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Observation of X(3872) in S —>■ X{ZS72)Ktz decays 
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We report the first observation of B” —^ X{3872){K'’'-k~) and evidence for X{3872){K‘’n^). 

We measure the product of branching fractions for the former to be B{B‘’ —>■ X{3872){K^ tt~ )) x 
B(X(3872) J/lpTT+n-) = (7.9±1.3(stat)±0.4(syst)) x 10"® and find that B° ^ X(3872)Js'* (892)° 
does not dominate the B° X(3872)_R'’'"7r~ decay mode. We also measure B{B^ —>■ 
X(3872)(ii'°7r+)) x B(X(3872) = (10.6 ± 3.0(stat) ± 0.9(syst)) x 10"°. This study 

is based on the full data sample of 711 fb“^ (772 x 10^BB pairs) collected at the T(45) resonance 
with the Belle detector at the KEKB collider. 

PACS numbers: 14.40.Pq, 12.39.Mk, 13.20.He 

About a decade ago, the Belle Collaboration discov¬ 
ered the A(3872) state [l| in the exclusive reconstruction 
of —1 A(3872)(^ J/'!/)7r+7r“)Ar+ Q. Considerable 
effort by both experimentalists and theorists has been 
invested to clarify its nature. As a resulC we know pre¬ 
cisely its mass (3871.69±0.17) MeV/c^ have a strin¬ 
gent limit on its width (less than 1.2 MeV at 90% confi¬ 
dence level) 0 and have a definitive assignment of 
1++ The A(3872) has been observed to decay to sev¬ 
eral other final states: J/ifu if'l 0) T/^7r+7r“7r° 0 


and 0, [^. The proximity of its mass to the 

]jO-f)*o threshold, along with its measured partial de¬ 
cay rates, suggests that it be a loosely bound “molecule” 
of 11° and D*’’ mesons 111 or an admixture of 
with a charmonium (cc) state [Hill. Some authors 
have advanced a QCD-tetraquark interpretation for the 
A(3872), and predict the existence of charged- and C- 
odd partner states that are nearby in mass |13l| . Exper¬ 
imental searches for charged- 0, [3] and C-odd [3 [3 
partners report negative results. However, since these 
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searches are restricted to states with narrow total widths, 
the published limits may not apply if the partner states 
access more decay channels and are thus broader. More 
experimental information on the production and decays 
of the X(3872) will shed additional light on its nature. 

In this paper, we present the results of searches for 
X(3872) production via the X(?,S12)K^'k~ and 

X(3872)A'g7r+ decay modes, where the X(3872) 
decays to The study is based on 711 fb“^ 

of data containing 772 x 10® BB events collected with 
the Belle detector 17| at the KEKB e+e” asymmetric- 
energy collider 18| operating at the T(45') resonance. 
In addition to selecting B —>■ X(3872)Ar7r signal events, 
the same selection criteria isolate a rather pure sample 
of i? —> ip'Kir events that are used for calibration. 

The Belle detector is a large solid-angle magnetic spec¬ 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero¬ 
gel threshold Cherenkov counters (ACC), a barrel-like ar¬ 
rangement of time-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter (ECL) composed of 
CsI(Tl) crystals. All these detector components are lo¬ 
cated inside a superconducting solenoid coil that provides 
a 1.5 T magnetic field. An iron flux return located out¬ 
side the coil is instrumented to detect AT® mesons and 
to identify muons (KLM). The detector is described in 
detail elsewhere [iTj. 

Monte Carlo (MC) samples are generated for each de¬ 
cay mode using EvtGen [13 and radiative effects are 
taken into account using the photos pa ckag e. The 
detector response is simulated using geant3 [21j . 

Charged tracks are required to originate from the in¬ 
teraction point (IP). To identify charged kaons and pi¬ 
ous, we use a likelihood ratio TZk/tt = + C,k), 

where the kaon (pion) likelihood Ck is calcu¬ 

lated using ACC, TOF and CDC measurements. For 
the prompt charged kaon (pion), we apply the criterion 
T^Kj-K {T^'njK) > 0.6. Here, the kaon (pion) identification 
efficiency is 93% (95%) while the probability of misiden- 
tifying a pion as a kaon (kaon as a pion) is 8% (7%). 
For the pion daughters from -tp' or A(3872), we impose 
T^-KjK > 0.4; the corresponding pion identification ef¬ 
ficiency is 99% and the misidentification probability is 
8%. Candidates for the ATg tt+tt” decay are formed 
from pairs of oppositely charged tracks having an invari¬ 
ant mass between 488 and 506 MeV/c^ (±4cr around the 
nominal mass of ATg). The ATg candidate is also required 
to satisfy the criteria described in Ref. [23 to ensure that 
its decay vertex is displaced from the IP. A track is identi¬ 
fied as a muon if its muon likelihood ratio is greater than 
0.1, where the muon and hadron likelihoods are calcu¬ 
lated by the track penetration depth and hit scatter in 
the muon detector (KLM). An electron track is identi¬ 
fied with an electron likelihood greater than 0.01, where 
the electron likelihood is calculated by combining dE/dx 
from the CDC, the ratio of the energy deposited in the 


ECL and the momentum measured by the SVD and the 
CDC, the shower shape in the ECL, ACC information 
and the position matching between the shower and the 
track. With the above selections, the muon (electron) 
identification efficiency is above 90% and the hadron fake 
rate is less than 4% (0.5%). 

We reconstruct J/tp mesons in the decay channel 

(£ = e or /r) and include bremsstrahlung photons that are 
within 50 mrad of either the or e~ tracks [hereinafter 
denoted as e“'"e“( 7 )]. The invariant mass of the J/ip can¬ 
didate is required to satisfy 3.00 GeV/c^ < < 

3.13 GeV/c2 or 3.06 GeV/c^ < M^+^- < 3.13 GeV/c^ 
(with a distinct lower value accounting for the residual 
bremsstrahlung in the electron mode). A mass- and 
vertex-constrained fit is performed to the selected J/'0 
candidate to improve its momentum resolution. The 
J/ip candidate is then combined with a tt+tt” pair to 
form an A(3872) [ip') candidate whose mass must sat¬ 
isfy 3.82 GeV/c^ < < 3.92 GeV/c^ (3.64 GeV/c^ 

< < 3.74 GeV/c^). The dipion mass must also 

satisfy — {mjj^ -I- 0.2 GeV/c^), where 

rrij/^ is nominal mass. This criterion corresponds to 
> 575 (389) MeV/c^ for the A:(3872) (^/>') mass 
region and reduces significantly the combinatorial back¬ 
ground [3 while also flattening the background shape dis¬ 
tribution in To suppress the background from 

e+e” —>■ qq (where q = u, d, s, c) continuum events, we 
require R 2 < 0.4, where R 2 is the ratio of the second- to 
zeroth-order Fox-Wolfram moments (^ . 

To reconstruct a neutral (charged) B meson can¬ 
didate, a AT+tt" [Kq-k^) candidate is combined with 
the A(3872) or ip' candidate. We select B candidates 
using two kinematic variables: the energy difference 
A A = Eb — Abeam and the beam-energy constrained 
mass Mbc = {y/E'^eam - 7 'bc 2 )/c^, where Abeam is the 
beam energy and Eb and pB are the energy and mag¬ 
nitude of momentum, respectively, of the candidate B 
meson, all calculated in the e+e” center-of-mass (CM) 
frame. Only B candidates having Mbc > 5.27 GeV/c^ 
and |AA| < 0.1 GeV are retained for further analysis. 
After all selection criteria, approximately 35% of events 
have multiple B candidates. For an event with more 
than one B candidate, we select the candidate having 
the smallest value of 


2 (Mbc - 5.2792 GeV/c2)2 , x| 
X — _2 + 


"Mbc 


ndf 


( 1 ) 


where UMbc is the Mbc resolution (estimated to be 
2.925 MeV/c^ from a fit to A® '>P'K^tt~ events), x% 
is the quality of the vertex fit of all charged tracks (ex¬ 
cluding the Kg daughters), ndf = (2A —3) in this fit and 
N is the number of fitted tracks. The correct candidate 
is selected in about 60% of the B A(3872)A7r events 
with multiple entries. 

To extract the signal yield of A —>■ A(3872)(—;• 
J/ip'K^'K~)K'K, we perform a two-dimensional (2D) un- 











4 


binned extended maximum likelihood fit to the /S.E and 
distributions. For the signal, the AE distribu¬ 
tion is parametrized by the sum of a Crystal Ball 241 
and a Gaussian function while the distribution 

is modeled using the sum of two Gaussians having a 
common mean. The 2D probability distribution func¬ 
tion (PDF) is a product of the individual one-dimensional 
PDFs, as no sizable correlation is found. 

The main background contribution in B 
(J/tjjTT~^TT~)K'K is expected to arise from inclusive B de¬ 
cays to J/V', which is confirmed by the low background 
found in the J/ijj mass sidebands in the data. To study 
this background, we use a large Monte Carlo sample of 
B —>• J/tpX events corresponding to 100 times the inte¬ 
grated luminosity of the data sample, and we find that 
few backgrounds are peaking in the distribu¬ 

tion (nonpeaking in the AE distribution) and vice versa. 
The remaining backgrounds are combinatorial in nature 
and are flat in both distributions. This background is 
parametrized by first-order Chebyshev polynomial. 

For the B^ X{3872)K'^tt~ decay mode, a 2D fit is 
performed. The mean and resolution of and AE 

are fixed for the X(3872) region from signal MC sam¬ 
ples after being rescaled from the results of the B^ 
ijj' decay mode. Further, we correct the mean of 
a Gaussian function for the distribution because 

of a difference in the shift of the ip' and X(3872) recon¬ 
structed and generated masses as seen in MG samples. 
The tail parameters, a and n of the Crystal Ball (CB) 
function, and the ratios of the CB’s area and width to the 
corresponding quantities of the Gaussian component are 
fixed according to the signal MC simulation. The peaking 
components can be divided into two categories: the one 
peaking in but nonpeaking in AE that comes 

from the B — X{3872)X' decays where the X(3872) de¬ 
cays in J/'ipTT^TT~ (here X' can be any particle), and the 
other peaking in AE but nonpeaking in which 

comes from a B with the same final state where 
is not from a X(3872). The peaking background in AE 
is found to have the same resolution as that 
of the signal, so the PDF is chosen to be the same as 
the signal PDF, while the nonpeaking background in the 
other dimension is parametrized with a first-order Cheby¬ 
shev polynomial. Parameters (slopes) of the background 
PDFs , which are of nonpeaking or combinatorial nature, 
are allowed to vary in the fit. The fits are validated on 
full simulated experiments and no significant bias is seen. 
Figured] (top) shows the signal-enhanced projection plots 
for the B° —>• X{3872){K^tt~) decay mode. The result 
of the fit and branching fractions derived are listed in Ta¬ 
ble [T| We find a clear signal for B^ —>■ X{8872)K^tt~ of 
116±19, signal events corresponding to a significance [2^ 
(including systematic uncertainties related to the signal 
yield as mentioned in TableHII of 7.0 standard deviations 
(cr), and measure the product of branching fractions to be 
S(B° ^ X(3872)iF+7r-) x S(X(3872) ^ J/V^Tr+Tr") = 


(7.9±1.3(stat)±0.4(syst)) x 10“®. The efficiency used for 
estimating the branching fraction is calculated from the 
individual efficiencies and the fractions of the different 
components obtained in the mass, as explained 

below. The statistical significance is estimated using the 
value of yJ—2 ln(£o/'^^max) where £max (^Cq) denotes the 
likelihood value when the yield is allowed to vary (fixed 
to zero). 

The above fit is validated on the calibration mode 
B^ —>■ Ip'. In contrast to the X(3872) mass region, 
the mean and resolution in both distributions 
and AE) are allowed to vary in the fit. Figure[T](bottom) 
shows the signal-enhanced projection plots for the 5° ^ 
Ip 'decay mode. We measure the branching frac¬ 
tion to be B{B^ ip'K'^TT~) = (5.79±0.14(stat)) x 10“^, 
consistent with an independent Belle result based on a 
Dalitz-plot analysis ( 2 ^ . 




FIG. 1: Projections of the {AE, fit for the ^ 

X(3872)(—>■ J/pTT'^'K~)'K~ decay mode (top) and the 
5° —> ip'{^ J'K~ decay mode (bottom): (a) AE 
distribution for 3.860 GeV/c^ < < 3.881 GeV/c^, (b) 

distribution for —11 MeV < AE < 8 MeV, (c) AE 
distribution for 3.675 GeV/c^ < < 3.695 GeV/c^ and 

(d) distribution for —11 MeV < AE < 8 MeV. The 

curves show the signal (red long-dashed curve) and the back¬ 
ground components (black dash-dotted line for the component 
peaking in but nonpeaking in AE, green dashed line 

for the one peaking in AE but nonpeaking in and 

magenta long dash-dotted line for combinatorial background) 
as well as the overall fit (blue solid curve). 


Further, to determine the contribution of the K*{892) 
and other intermediate states, we perform a 2D fit to 
AE and in each bin of Mk^ (0.1 GeV/c^ wide 

bins of Mktt in the range [0.62,1.42] GeV/c^), which 
provides a background-subtracted Mktt signal distribu- 
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TABLE I: Signal yield (Y) from the fit, weighted efficiency (e) after particle-identification correction, significance (E) and 
measured B for X{3872)K+'k- and B+ X{3872)K°n+. The first (second) uncertainty represents a statistical 

(systematic) contribution. 


Decay mode 

Y 

e{%) 

S (cr) B{B 

^ A(3872)A'7r)x S(A(3872) -)• Jlip-K+ir-) 

5° ^ A(3872)A:+7r“ 

116 ± 19 

15.99 

7.0 

(7.9 ±1.3 ±0.4) X 10“° 

B+ ^ A(3872)A:°7r+ 

35 ±10 

10.31 

3.7 

(10.6 ±3.0 ±0.9) X 10“® 


TABLE II: Summary of the systematic uncertainties in per¬ 
cent. 


Source 

A(3872) A(3872) 


A:+7r“ A:°7r+ 


Lepton ID 

3.4 

3.4 

Kaon ID 

1.1 


Pion ID 

2.5 

3.2 

PDF modeling 

+ 1.8 
-1.3 

+ 4.2 

- 4.9 

Tracking efficiency 

2.1 

2.5 

K% reconstruction 


0.7 

^BB 

1.4 

1.4 

Secondary B 

0.4 

0.4 

Efficiency 

0.6 

1.0 

Fit bias 

0.6 

3.1 

Total 

5.4 

8.0 


tion. All parameters of the signal PDFs for 
and AE distributions are fixed from the previous 2D 
fit to all events. We perform a fit to the Mktv dis¬ 
tribution using AT* (892)° and (Ar+7r“)NR components, 
which are histogram PDFs obtained from MC samples. 
Note that the ^ A(3872)Ar2*(1430)° decay is kine¬ 
matically suppressed. We do not consider the inter¬ 
ference between the AT* (892) and nonresonant compo¬ 
nent since the number of candidates is not large enough 
to make a full amplitude analysis. The resulting fit 
result is shown in Figure [UJa). We obtain 38 ± 14 
(82 ± 21) signal events for the —>■ X(3872)Ar*(892)° 

(5° A(3872)(A'+7r- )nr) decay mode, whose sum is 

consistent with the total yield obtained from the global 
fit. This corresponds to a 3 .Oct significance (includ¬ 
ing systematic uncertainties related to the signal yield) 
for the f A'(3872)(—>• J/'ipTr^Tr~)K*{892)^ decay 
mode, and a product of branching fractions of B{B° —>• 
A(3872)A:*(892)0) x S(A(3872) JltjjTT+T:-) = (4.0 ± 
1.5(stat) ±0.3(syst)) x 10“°. The ratio of branching frac¬ 
tions is 

B{B° X(3872)A:*(892)0) x B{K*{S92f ^ AT+Tr”) 
B{B^ X(3872)A:+7r-) 

= 0.34 ± 0.09(stat) ± 0.02(syst). 

( 2 ) 

In the above ratio, all systematic uncertainties cancel ex¬ 



FIG. 2: Fit to the background-subtracted Mktt distribution: 
(a) for the B° —>■ A (3872) (A'”''tt”) decay mode, the curves 
show the B^ —^ A(3872)A'*(892)° (red long-dashed lines), 
B° A(3872)(A'+7r“ )nr (green dot-dashed lines), as well as 
the overall fit (blue solid lines), (b) for the ijj' 

decay mode, the curves show the 5° —>■ ip'K* (892)° (red long- 
dashed lines), B° —>■ ip '(green dot-dashed lines), 
^ Ip'K 2 (1430)° (magenta dashed lines) as well as the 
overall fit (blue solid lines). 


cept those from the PDF model, fit bias and efficiency 
variation over the Dalitz distribution. We generate pseu¬ 
doexperiments to estimate the significance of the fit. 

The B^ ip' mode is analyzed with the same 
procedure, with 0.051 GeV/c^ wide bins, due to the copi¬ 
ous yield, and in the Mkit range [0.600, 1.569] GeV/c^. 
We perform a fit to the obtained Mk-k signal dis¬ 
tribution again to extract the contributions of the Ktt 
nonresonant and resonant components. For this pur¬ 
pose, we use histogram PDFs obtained from MC sam¬ 
ples of several possible components of the (AT+tt") sys¬ 
tem: Ar*(892)°, Ar2(1430)° and nonresonant AT+tt” 
((Ar+7r“)NR); in the last case, 5° '0'(Ar+7r“)NR is gen¬ 

erated uniformly in phase space. The fit result is shown 
in Figure [2l[b). The AT* (892) dominates and we measure 
B{B° -s- iP'K* (892)°) = (5.88±0.18(stat)) x 10“^, which 
is consistent with the world average Q. 

In contrast to ip'{K'^7r~) [where the ra¬ 
tio of branching fractions is 0.68 ± 0.01 (stat)], —>• 

A(3872)Ar*(892)° does not dominate in the 5° 
A(3872)A:+7r-. 

We also investigate the decays 5+ ^ A(3872)(^- 
J/'(/i7r+7r“)(Ar°7r+). The PDFs of AE and are 

the same as those for the neutral B meson decay mode. 
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The projections of the 2D fit for —>■ X(3872)(^> 

J/'(/i7r+7r“)(i^°7r+) in the signal-enhanced regions are 
shown in Figs. EJa) and EJb). We find 35 ± 10 events for 
the X(3872)(—>• J/'07r+7r“)(if°7r+) decay mode, 

corresponding to a 3.7 (t significance (including system¬ 
atic uncertainties). The product of branching fractions is 
B{B+ X{m2)K^TT+) X S(X(3872) ^ = 

(10.6 ± 3.0(stat) ± 0.9(syst)) x 10“®. The above fit is 
validated for the i/j' mass region. The projections of the 
2D fit for B'^ —>• J/ilnr~^TT~)(K^TT~^) in the signal- 

enhanced regions are shown in Figs. [3Kc) and [3Kd). The 
branching fraction for —>■ J/'07r+7r“)(iF°7r+) is 

(6.00±0.28(stat)) x 10“'^, while the world average of this 
quantity is (5.88 ± 0.34) x 10“"^. 






FIG. 3: Projections of the {AE, fit for the —>■ 

X(3872)(—>■ J/ip-K'^'K~)Kg-K^ decay mode (top) and for the 
—>■ Ip'J/'tpTT'^'K ~)decay mode (bottom): (a) AE 
distribution for 3.859 GeV/c^ < < 3.882 GeV/c^, (b) 

distribution for —11 MeV < AE < 9 MeV, (c) AE 
distribution for 3.675 GeV/c^ < < 3.695 GeV/c^ and 

(d) distribution for —11 MeV < AE < 9 MeV. The 

curves show the signal (red long-dashed cnrves) and the back¬ 
ground components (black dash-dotted curves for the compo¬ 
nent peaking in bnt nonpeaking in AE, green dashed 

lines for the one peaking in AE but nonpeaking in 
and magenta long dash-dotted lines for combinatorial back¬ 
ground) as well as the overall fit (bine solid lines). 


Equal production of neutral and charged B meson 
pairs in the T(4iS') decay is assumed. Secondary branch¬ 
ing fractions used for calculation of B are taken from 
Ref. Q . Systematic uncertainties are summarized in Ta¬ 
bles El and mil A correction for small differences in the 
signal detection efficiency between signal MC events and 
data due to lepton, kaon and pion identification differ¬ 
ences is applied; samples of J/ip and ^ 


TABLE III: Summary of systematic uncertainties (in per¬ 
cent) used for the Mktt background-subtracted fit in B° —>■ 
A(3872)A+7r-. 


Source 

A(3872)A*(892)“ 

Lepton ID 

3.4 

Kaon ID 

1.1 

Pion ID 

2.6 

PDF modeling 

+ 1.5 
- 1.4 

Tracking efficiency 

2.1 

Kbb 

1.4 

Secondary B 

0.4 

MC statistics 

0.2 

Fit bias 

4.6 

Total 

7.0 


D°(—>■ Ar“7r+)7r+ decays are used to estimate this correc¬ 
tion. The uncertainties on these corrections are included 
in the systematic error. The uncertainty due to the fitting 
model is obtained by varying all fixed parameters by iltr 
and considering the corresponding change in the yield as 
the systematic error. The uncertainties due to tracking 
efficiency. Kg reconstruction and are estimated to 
be 0.35% per track, 0.7% and 1.4%, respectively. The 
systematic uncertainty due to secondary branching frac¬ 
tions is included. The systematic uncertainty for the ef¬ 
ficiency arises from the limited MC statistics and the ef¬ 
ficiency variation over the Dalitz distribution is also con¬ 
sidered. Small biases in the fitting procedure, estimated 
in the ensemble study, are also considered as a source of 
systematic uncertainty. For this study we perform a fit 
to 100 pseudoexperiments after embedding signal events 
obtained from MC samples to each inclusive MC sam¬ 
ple. All the above stated systematic uncertainties are 
added in quadrature and result in a total systematic un¬ 
certainty of 5.4%, 8.0%, 7.0% for B^ X(3872)A:+7r-, 
B+ A(3872)A:^7r+ and A(3872)A:*(892)°, re¬ 

spectively. 

In summary, we report the first observation of 
the A(3872) in the decay -)> X(3872)A:+7r-, 
A(3872) ^ J/V'tt+tt-. The result for the A(3872), 
where B^ —)> A(3872)A*(892)° does not dominate 
the B^ —>• A(3872)(A+7r“) decay, is in marked 

contrast to the ip' case. We have checked for a 
structure in the A(3872)7r and A(3872)A invariant 
masses and found no evident peaks. We mea¬ 
sure S(B° A(3872)(A+7r-)) x S(A(3872) ^ 

J/ipTT^TT~) = (7.9 ± 1.3(stat) ± 0.4(syst)) x 10“® and 
B{B+ A(3872)A°7r+) x S(A(3872) ^ JIiPtt+tt-) = 
(10.6 ± 3.0(stat) ± 0.9(syst)) x lO"®. 
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